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ABSTRACT 

We present the results of a deep, wide-field search for transiting 'Hot Jupiter (HJ)' planets in the 
globular cluster lo Centauri. As a result of a 25-night observing run with the ANU 40-inch telescope 
at Siding Spring Observatory, a total of 109,726 stellar time series composed of 787 independent data 
points were produced with differential photometry in a 52'x52' (0.75 deg^) field centered on the cluster 
core, but extending well beyond. Taking into account the size of transit signals as a function of stellar 
radius, 45,406 stars have suitable photometric accuracy (<0.045 mag to V=19.5) to search for transits. 
Of this sample, 31,000 stars are expected to be main sequence cluster members. All stars, both cluster 
and foreground, were subjected to a rigorous search for transit signatures; none were found. Extensive 
Monte Carlo simulations based on our actual data set allows us to determine the sensitivity of our survey 
to planets with radii '^1.5Rjup, and thus place statistical upper limits on their occurrence frequency F . 
Smaller planets are undetectable in our data. At 95% confidence, the frequency of Very Hot Jupiters 
(VHJs) with periods P satisfying Id < P < 3d can be no more than i^vHj < 1/1040 in uj Cen. For HJ 
and VHJ distributed uniformly over the orbital period range Id < P <5d, Fvhj+h,7 < 1/600. Our limits 
on large, short-period planets are comparable to those recently reported for other Galactic fields, despite 
being derived with less telescope time. 



Subject headings: globular clusters: individual (NGC 5139, lo Centauri) 
techniques: photometric 



planetary systems 



1. INTRODUCTION 

The identification and subsequent study of extrasolar 
planets has become a subject of intense interest in recent 
years. To date, '^250 giant planets have been discoverecfl, 
ranging in mass from Neptune to greater than Jupiter. 
These new worlds are altering our understanding of the 
formation and evolutionary processes of giant planets in 
the immediate solar neighborhood. Currently, the major- 
ity of them have been found using radial velocity (RV) 
techniques, which favours the detection of close-in, mas- 
sive planets. 

Statistical analysis of current RV detections indicates 
that 1.2%±0.3% of nearby F,G and K stars are orbited by 
'Hot Jupiter' (HJ) planets, those with orbital periods of 
only a few days (<0.1 AU) and minimum masses a pproxi- 
mately equal to that of Jupiter ( Marcv et al.ll2005h . Such 
discoveries have challenged traditional ideas of planetary 
evolution, implying that a rapid migration of the planet 
takes place soon after formation. 

Planet frequency from radial velocities appe ar to de 
pend on the metallicity of the host star (iGonzale: 
igpT*; 'LaughlinlllOOS ISantos et al.][200l ': 'Fischer & Valent 
2005). However, there is very little observational evidence 
for a lower planetary frequency at quite low metallicities, 
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due to the bias of radial velocity detections in association 
with bright nearby high metallicity stars. Hence dedicated 
surveys in low metallicity environments, like globular clus- 
ters, can provide information to help understand this re- 
lationship in a more robust manner. Indeed, does low 
metallicity halt planet formation or just affect the plane- 
tary migration process? 

Of the currently known planets, ^^90% have only a min- 
imum mass assigned to them, due to the unknown incli- 
nation of the planetary orbit, and unknown radii and den- 
sities. These quantities can be measured for transiting 
planets. Due to its short orbital period, each HJ has a 
non-negligible probability of transiting its host star that 
depends on the orbital separation and the ratio between 
the stellar and planetary radii. Typical transit depths and 
durations are ~1.5% and ^2 hours respectively. With pre- 
cise photometry these transiting systems can be identified 
in the field, leading to direct measurements of the plan- 
etary radius (provided the stellar radius is known) from 
the depth of the transit dip. Studies of the planetary at- 
mosphere may be attempted if transits occur, and when 
coupled with RV measurements, accurate mass and den- 
sity determinations can be made. 

Currently, 28 tran siting ex oplanets ar e known 
(|Charbonneau et al.l l200ol: iHenrv et al.] 120001: 



^ http://exoplanet.eu/ 
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discovered through radial velocity searches. Despite an in- 
creasing number, the detection rate of planets from transit 
search e s is s ignificantly lower than initially expected eg, 
iHornd ()2003[ ). This lack of detections is due, in part, to 
observing strategy: a long observing window (>1 month 
equivalent) with a dedicated telescope coupled with a wide 
field and high temporal resolution are needed to sample 
enough stars frequently enough to allow the detection of 
a transit. Also the production of a large enough number 
10^ -> 10^) of high-quahty (<0.02 mag) lightcurves of 
dwarf stars coupled with a sensitive detection algorithm 
with low fal se alarm rates is no n-trivial. 

Recently, iGould et al.] (|2006D analyzed OGLE III tran- 
sit surveys in Galactic fields and concluded that the occur- 
rence frequencies of the detected planets in these surveys 
is not statistically different from that found in RV sur- 
veys of nearby stars. However, they did conclude that 
the frequency of HJ planets with periods P satisfying 
3<P<5 days {F = l/320(llj;g) at a 90% confidence) 
was statistically different from and larger than that of 
Very Hot Jupiter (VHJ) planets (1<P<3 days), which 
have F =l/710(llo5!^) at 90% confidence. Since the 
OGLE HI surveys detected no planets with radius larger 
than 1.3Rjup,they placed upper limits on the occurrence 
frequency of larger worlds. 

Transit searches, unlike RV, are not limited to the im- 
mediate solar neighborhood and can be used to measure 
relative planet frequencies in various regions of the Galaxy, 
providing information on the role environmental effects 
play on HJ planet formation. Early predictions for the suc- 
cess o f trans it searches in open clusters was presented by 
iJanesI (|1996l ). indicating that with a large amount of tele- 
scope time planets could be detected via the transit tech- 
nique in nearby clusters. More recently, iPepper fc Gaudil 
(|2005[ ) presented an analysis of the prospective harvest 
of cluster transit surveys by discussing the observational 
techniques and methods to maximise the chances of a 
detection. They concluded that due to their mass func- 
tions, the most populous, nearby and bri ght clusters have 
the greatest chance of yielding a planet. iPepper fc Gaudil 
((2OO6I then went on to discuss specifically the detection of 
short period 'Hot Earth' and 'Hot Neptune' planets. The 
detection yields for various nearby clusters with v arious 
instruments was estimated. lAigrain fc PontI (|2007D agree 
that small- aperture wide field surveys targetting nearby 
clusters have the potential to discover transiting Hot Nep- 
tune planets. 

Transit searches have been undertake n in bright 
meta l-rich open cl usters including STEPPS (iBurke et al.l 
[2006'). USt AP S (IStreet et a l.' '2003': 'Hood ct al.' '2005), 
PISCES (.Mocheiska et all 12003. 2005 . 2006). Monitor 
(lAigrain et al.l 120071). EXPLORE-OC (jvon Braun et al.l 
I2OO5D and liT the Praesepe cluster {M44) with KELT 
( Pepper et al.ll200"7t ). Searches have also been per f ormed 
in the general Galactic field (jUdalski et al.l l2002l 120041: 



iHidas et al.ll2005l : iKane et al.ll2005l: IWeldrake et alll2007D 
and toward the Galactic Bulge (jSahu et al. I I2006D . If clus- 
ter candidates are confirmed as planetary in nature, dif- 
ficult if fainter than V^17.0, they can provide informa- 
tion on the timescales of HJ formation and subsequent 
migration. Null results of high significance allow planet 
frequency upper limits to be estimated. 

Globular clusters provide an excellent opportunity to 
study the effects of environment on planetary frequency. 
Two bright, nearby southern clusters, 47 Tuc and lu Cen- 
tauri, have stars in sufficient numbers (~10^) and bright- 
ness (V^17) for meaningful statistics to be gained us- 
ing ground-based telescopes of moderate aperture. 47 
Tuc was previously sampled for planetary transits, result- 
ing in a hig h significance (>3a) null result in both the 
cluster core ("Gilliland et al.l (20001 ) and in the outer halo 
( Weldrake et al. 2005). These two results strongly indicate 
that system metallicity - not crowding - is the dominant 
factor determining HJ frequencies in this cluster. 

This paper presents the results of a dedicated transit 
search in the second cluster, uj Centauri, in order to test 
further the dependence of planetary frequency on stellar 
metallicity and crowding. Omega Centauri has only 1/10**^ 
the core density of 47 Tuc yet contai i is five times the to- 
tal mass (5.1x10*^Mq, iMevlan et all (|l995t) ). Due to its 
low stellar density compared to other globular clusters and 
long stellar interaction timescale, a null result for uj Cen- 
tauri can be used to test the relative importance of stellar 
metallicity over density in the formation of giant planets. 

Omega Centauri {oj Gen, NGC 5139) has been subjected 
to intense research over the years. The cluster is unique 
among globular clusters in that i t displays a distinct spread 
of metallicity among its stars ([Dickens fc Woollev 1967 



Norris fc Bessen"1975t lLee et al.il999l:lPancino et al.lboOO : 

SpUinia et al. 2005), due to an extended period of star for- 
mation and chemical enrichment. Using He abundances, 
iNorrisI (120041 ) has shown that the cluster has three dis- 
tinct stellar populations, with metallicities of —1.7, —1.2 
and —0.6 dex, corresponding to 0.80, 0.15 and 0.05 of the 
total population respectively. 

The cluster has a highly-bound, retrograde orbit 
(jDinescu et al.lll999D and is bv far the most massive of the 
globular clusters (jMevlan et al.lll995D . Indeed, these va- 
garies have led to the theory that the cluster had an exter- 
nal origin, bein g the left-over rema i ns of a tidally disrupted 
dwarf galaxv (iBekki fc FreemanI 120031 : lldeta fc Makinol 
12004 iBekki fc Norrisll2006D . With its relative proximity, 
u> Gen presents a statistically significant number of upper 
main sequence stars that can be searched for transiting HJ 
planets. 

Here we present the result from a vigorous search for the 
transit signatures of large planets on 45,406 lightcurves in 
a 0.75deg^ field centered on uj Gen. The same set of ob- 
servations yielded a total of 187 variable stars in the field, 
81 of which are n ew discoveries, and are presented in a 
companion paper ([Weldrake et al.|[2007bD . Furthermore, 
we observed a control field in the Lupus Galactic Plane to 
test the data reduction and transit identification strate- 
gies. Analysis for this field is ongoing, but has led to the 
identification of several transit candidates, of which none 
similar were seen in the to Gen dataset. One candidate 
in particular has excellent prospects for being a new Hot 
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Jupiter planet (jWeldrake et al.ll2007l ). 

Section 2 of this paper describes our observational strat- 
egy and data reduction details. Section 3 details how the 
photometry was obtained for both the crowded core re- 
gions and outer halo parts of the dataset. The cluster 
Color Magnitude Diagram dataset (along with astrome- 
try) are also briefly discussed. Section 4 describes the 
stellar parameters of the cluster stars that were searched 
for transits, and the expected characteristics of the tran- 
sits themselves. The total number of stars in the field 
(both in the cluster and the foreground galactic disk) is 
also calculated. Section 5 describes the transit detection 
algorithms used in our search and our removal of system- 
atics in the photometry. Our Monte Carlo simulations to 
derive expected recovery of real transits and false alarms 
are described in Section 6, and their application to esti- 
mate our HJ sensitivity outlined in Section 7. The results 
of our transit search in w Cen are presented in Section 8, 
with discussion, comparison to the literature and interpre- 
tation in Section 9. We conclude in Section 10. 

2. OBSERVATIONS AND DATA REDUCTION 

Our dataset was produced using the Australian National 
University (ANU) 40-inch (Im) telescope located at Sid- 
ing Spring Observatory, fitted with the ANU Wide Field 
Imager (WFI). This telescope and detector combination 
permits a 52'x52' (0.75 deg^) field of view, which was 
centered on the cluster core. 

WFI is capable of sampling a large fraction of the clus- 
ter in a single exposure, as the field extends to 50% the 
cluster tidal radius (Harris 1996), allowing a large time 
series dataset to be produced with only a single pointing. 
The WFI detector consists of a 4x2 array of 2048x4096 
pixel back-illuminated CCDs in an 8Kx8K arrangement. 
The pixel scale at the telescope Cassegrain focus is 0".38 
pixel" ^, permitting suitable sampling of the point spread 
function (PSF). 

In order to increase the likelihood of transit detection, 
an observing strategy was employed to maximize the tem- 
poral resolution of the dataset while keeping the resultant 
signal-to-noise ratio (S/N), and hence photometric preci- 
sion, high. In other words, deep images with short expo- 
sure times are necessary. In order to do this with a Im 
telescope, a special broad-band filter was constructed to 
cover the combined wavelength range of Cousins V and R. 

A five minute exposure of a V=18.5 star (typical of a 
main sequence star in the cluster) with this filter has a 
photon noise S/N of 220 in 7-day moon and 2" seeing. 
In order for a typical HJ (transit depth ^^0.015 mag) to 
be detected at the 3a level, an effective S/N of 200 is re- 
quired, which is thus obtainable with this V-l-R filter on 
timescales that would well resolve the expected '^2-hour 
transits. This same telescope and detector combination 
was used successfully for a transit search in 47 Tucanae, 
which led to a high si gnificance null result in that cluster 
(IWeldrake et al.ll2005D . 

The globular cluster lu Centauri was observed for 25 
contiguous nights, from 2003 May 2 to 2003 May 27 
with a field center (J2000.0) of a = 13''26'"45.89^ S = 
-47°28'36.7" (J2000.0). An exposure time of 300s was 

^ IRAF is distributed by National Optical Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc., 
under cooperative agreement with the National Science Foundation. 



used throughout the observing run with the V-l-R filter. 
Each image was checked individually for quality after read- 
out at the telescope. If an image displayed bad seeing 
(>3"), bad focus or cosmetic problems (such as satellite 
trails, intermittent clouds or other adverse effects) it was 
discarded. After this quality control, a total of 875 images 
were obtained over the 25-night run, the dataset having 
an average temporal resolution of six minutes and cover- 
ing nine hours on a typical night. 

Initial image reduction was performed using standard 
reduction practices with the MSCRED package of IRAfQ. 
This procedure included region trimming, overscan cor- 
rection, bias correction, flat-fielding and dark current sub- 
traction. The reduced images were then checked for qual- 
ity before entering the photometric pipeline. Of the 875 
images obtained, 90% (787 images) were deemed suitable 
for use in the production of the time series dataset, as 
indicated by their small telescope offsets and good seeing. 

3. PHOTOMETRY AND PHOTOMETRIC ACCURACY 

High precision photometry can be obtained on faint tar- 
gets in the crowded field of lo Cen by performing differ- 
ential photometry. This method was originally described 
as an opti mal Point-Spread-Funct ion (PSF) matching al- 
gorithm bylAlard & Lupton] (|1998D . and was subsequently 
modified by Wozniak (2000) for use in detecting microlens- 
ing events. A detailed description of the Difference Image 
Analysis (DIA) method and software pipeline can be found 
in Wozniak's paper; only the main steps are summarized 
here. 

The process of matching the stellar PSF in a large 
database of images reduces dramatically the systematic ef- 
fects due to varying atmospheric conditions on the photo- 
metric precision, allowing the detection of small brightness 
variations in faint targets with ground-based observations. 
DIA is also an excellent photometric method for dealing 
with crowded fields. Since a larger number of pixels con- 
tain information on any PSF differences as the number 
of stars increases, there is an improvement in the PSF- 
matching process. Baseline fiux measurements of the stars 
are made via profile fitting on a master template frame, 
which is produced by median-combining a large number 
(40-1- in our case) of the best quality images with small 
offsets. This fiux measurement is used as the zero-point in 
the photometric time series for an individual star. 

Stellar positions are determined from a reference image, 
usually the image with the best seeing conditions; and all 
other images in the dataset (including the template) are 
shifted to match. The best PSF-matching kernel is then 
found for each image, and each registered image is sub- 
tracted from the template. The residual subtraction im- 
age is generally dominated by photon noise. Any object 
that has changed in brightness between the image and the 
template is given away as a bright or dark spot. 

The pixel coordinates of all visible stars were determined 
separately on the reference frame via DAOFIND within 
IRAF, and the profile photometry was then extracted from 
the subtracted frames at those determined positions. 

Differential photometry produces a time series measured 
in differential counts, a linear flux unit from which a con- 
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stant reference flux (taken from the template) has been 
subtracted. To convert this to a standard magnitude sys- 
tem, the total number of counts for each star was measured 
using the PSF photometry package of DAOPHOT within 
IRAF, with the same images and parameters as used in 
the photometry code. The time series fluxes were then 
converted into magnitude units via the relation: 

Am, = -2.51og[(7V, + 7V,.cf,^)/iV,ef,.] 

where iVrof.i is the total flux of star i on the template 
image and Ni is the difference flux in the time series as 
produced with the photometric code. 

When combining differential fluxes with DAOPHOT- 
derived photometry on the reference image, it is impor- 
tant to correct for errors based on the individual apertures 
used. The scaling between the two fluxes was determined 
via performing an aperture correction on the DAOPHOT 
magnitudes for the cluster stars, as per the method de- 
scribed in Appendix B of Hartman ct al. (2004). Wc found 
that our PSF magnitudes were consistently 0.09 magni- 
tudes brighter than the aperture-derived values (using the 
same aperture values as in the differential photometry). 
We corrected for this by shifting our magnitude zero-point 
to 25.0- 0.09 = 24.91. 

A total of 109,726 stellar photometric time series were 
produced across the whole WFI fleld, each containing 787 
independent data points, which then became the subject 
of analysis. 

3.1. Photometry of the cluster core 

Data from each of the four outer CCDs of WFI were 
divided in half for DIA analysis, with each half producing 
an average of 9,500 time series. For the crowded core of 
the cluster, where the number of stars becomes very large, 
computational limitations necessitated a different strategy. 
For the core regions, the images were analysed in 360 in- 
dividual subframes, 90 per CCD. The locations of these 
subframes were chosen so that no stars would be lost at 
the edges of each individual subframe, and the entire core 
region of the cluster was covered. Data of sufficient qual- 
ity could not be obtained in regions affected by telescope 
offsets during observing (typically a 160 pixel border sur- 
rounding each C CD. The SYSREM systematics removal 
package of lTamuz et al.l (|2005[ ) was applied to all resulting 
lightcurves. 

Fig. [T] presents the resultant DIA-derived photometric 
precision, measured as root-mean-square scatter (rms), for 
the 97,935 stars that were cross-identified with the cluster 
CMD dataset. The position of the cluster main sequence 
turnoff (MSTO) is marked to indicate where the cluster 
stars become members of the main sequence. Objects to 
the left of this line are likely Red Giant Branch and fore- 
ground Galactic Disk stars. The photometric uncertainty 
is 1% at V=17.4, rising to 4% at V=19.0. 

Also overlaid on Fig. [l]is the expected depth of a tran- 
sit for a 1.5Rjup planet as a function of V magnitude of 
the parent star. The total star, background and residual 
noise contribution is also plotted (thin locus of red points) . 
The photometry is photon-noise-limited for V^17.0 and 
sky-t-residuals dominated to fainter magnitudes. We de- 
fine our lower limit as lying at the magnitude value where 



all of the stars in our dataset have photometric uncertainty 
larger than the expected transit depth. It can be seen on 
Fig. [l]that this limit is reached at V=19.5. 

3.2. Color Magnitude Diagram and Astrometry 

Using the same telescope/detector combination, a V, V- 
I color magnitude diagram (CMD) totalling 203,892 stars 
was produced for the observed field, enabling the detected 
transiting systems to be placed on the standard V and I 
magnitude system. This aids in determining their likely 
nature. The total CMD dat aset is presented in ou r vari- 
able star companion paper (iWeldrake et al.ll2007bn . Fig. 
[2] presents the diagram for the upper main sequence region 
that was targeted for the transit search. The photomet- 
ric errors reported by DAOPHOT in both V and V-I are 
marked as errorbars as a function of V magnitude. 

As standard field data was of unacceptable quality, the 
CMD was calibrated b y matchin g stella r astrometry from 
our catalog to that of IColemanI (|2004| ) (also taken with 
the ANU 40-inch and WFI combination in V and I). The 
diff erence in V and I between our uncalibrated data and 
the lColemanI (j2004[ ) calibrated data was measured for each 
of the matched stars (totalling more than 20,000) in each 
CCD independently. The resulting calibration accuracy is 
<0.03 magnitudes. Variable stars found in the monitoring 
data were identified on the CMD by visually identifying 
the star on the template image and comparing it to the 
V-band CMD image for the respective CCD. 

Astrometry was obtained for a total of 212,959 stars 
identified in the V band image of the cluster and 243,466 
stars in the I-band. A search of the USNO CCD As- 
trograph Catalog (UCACl) was carried out for astromet- 
ric standard stars within the field. Several hundred such 
stars were successfully identified, producing an accurate 
determination of the astrometric solution for the stars in 
each CCD independently; the resulting calibration accu- 
racy was 0.25". The extent of our single WFI field of 
view can be seen in Fig. [31 plotted as ARA and ADec in 
degrees from the position of the cluster core. For compari- 
son, the locations of the cluster core radius (inner ellipse), 
half-mass radius (middle ellipse), and the location of 50% 
of the cluster tidal radius are also marked overlying the 
location of the eight WFI CCDs. 

4. STELLAR PARAMETERS, TRANSIT EXPECTATIONS AND 
TOTAL NUMBER OF STARS 

In order to determine the expected depths and dura- 
tions of planetary transits (which are necessary for the 
Monte Carlo simulations) against main sequence cluster 
stars, the radii of those stars were first calculated as a 
function of magnitude. Th is was done by producing three 
theoretical I Yi et al] (j2003l ) isochrones, each with metallic- 
ity and age values for uj Cen as taken from Norri^ (120041 ) 
to simulate the cluster stellar content. The majority of the 
stars (~80%) are expected to lie on the most metal-poor 
([Fe/H]=— 1.7) sequence, which was used for our analysis. 
These three isochrones and the corresponding parameters 
can be seen superimposed on Fig. [2l 

Table 1 presents the resulting stellar parameters as de- 
termined from the isochrone for this metal poor popula- 
tion for stars with V=17. 0^20.0, namely those stars on 
the cluster turnoff and upper main sequence. Tabulated 
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are the apparent magnitude, the stellar mass, temperature 
and luminosity (in Solar units), the stellar surface gravity 
and the corresponding absolute magnitude. The stellar 
radius has been calculated from the surface gravity values 
and is also tabulated. The stellar radius increases rapidly 
for cluster stars brighter than the main sequence turnoff 
(at V^17.5), hence transit visibility becomes reduced for 
the brightest stars. 

From these calculated radii, the expected depths in mag- 
nitude units (Dep) and durations in hours (Dur) for a tran- 
siting planet were found via: /„ ^ 2 



4.1. Galactic Disk Contamination 



(1) 



where Rp is the radius of the planet and is the radius 
of the star, and: 

Dur 1.412M::^/^R»P;^/^ (2) 

where and i?, are the stellar mass and radius in solar 
units, and Pnrh is the orbital per iod of the planet in days, 
taken from lGilliland et aD (|2000( ). We assume a period of 
3.3d, typical of HJ planets found in the solar neighbor- 
hood. The duration value incorporates a 7r/4 reduction 
in transit duration, corresponding to the average chord 
length across a stellar disk. A centrally-crossing transit 
would have a duration 1.28 times longer than the values 
determined here. This relation can be used to infer the 
planetary radius of any detected candidates with a known 
period and total duration. The orbital inclination is also 
determined in that calculation. 

The transit depth and duration values as calculated via 
the Yi isochrones are displayed for planets of l.SRjup in 
Fig. [5] for various V magnitudes of the cluster main se- 
quence. Also overplotted is the location in brightness 
at which various photometric precisions can be expected. 
The precision becomes comparable to the expected depth 
of a 1.5Rjup planet at V'^IQ.S. 

The expected transit depths are superimposed on Fig. 
[T]for all cluster stars. The hue indicates l.SRjup planet 
depths and we do not have the photometric accuracy to 
search for planets smaller than this limit. The transit 
depth becomes very small once the stars leave the clus- 
ter turnoff (hence becoming physically larger) and rapidly 
become undetectable. For larger planets, the transit depth 
is similar to the photometric uncertainty, as stars are 
smaller further down the main sequence. The optimal zone 
for transit detection is V=17.0— >19.0, with the limit at 
V~19.5, where photometric errors dominate the lightcurve 
and the transit depth is greater than the photometric un- 
certainty. We do not consider the signal-to-noise of the 
transits in our limit determination, only the magnitude 
limit for which the per-data point photometric error be- 
comes larger than the expected transit depth. 

With this knowledge, the total number of stars avail- 
able for analysis was determined for this planetary radius 
limit. For l.SRjup planets, 45,406 appropriate stars are 
in our dataset between cluster turnoff and V=19.5. How- 
ever, this number also contains some Galactic contamina- 
tion. From applica tion of a Bccanson Galactic field model 
(|Robin et al.l l2003[ ) for a WFI field in the direction of uj 
Gen, this contamination has been estimated at 31%. By 
accounting for Galactic contamination in this magnitude 
range, we arrive at ■~31,000 cluster stars suitable for a 
search for large transiting planets. 



Using the Becanson model ()Robin et al.ll2003l ). we can 
estimate the total number of foreground Galactic disk stars 
in the field (non cluster members) that have photometric 
uncertainty <2%. Assuming these stars are foreground 
stars with solar radius or lower, they are suitable for a 
transit search in their own right. Gluster stars of this 
brightness correspond to subgiant and red giant branch 
members, which are unsuitable for the search. 

There are an estimated total of 6500 foreground disk 
stars in the field as determined from the model with <2% 
photometry to a magnitude limit of V=18.2. Of these, 
770 are G type main sequence stars with 160 K and 25 M 
main sequence stars. The remainder are of unsuitable lu- 
minosity class. By considering the small radius of M-dwarf 
stars, it is possible to extend a search to V=19.5 with 
sufficient photometric uncertainty to detect a planet with 
radius approximately equal to Jupiter. This permits a to- 
tal estimate of ^^200 foreground M-dwarfs in the dataset. 
Although these numbers are insufhcient for a statistical 
transit detection, all stars were nevertheless searched in 
the course of this work. 



4.2. Dataset Systematics Removal 

Before being subjected to the main transit search, the 
time series dataset was filtered for outlying data points 
and subjected to the systematics removal package of 
iTamuz et al.l ((lOO^). Systematics are well known as the 
source of the vast majority of false transit-like features 
detected in large datasets, caused by various uncontrol- 
lable factors inherent in ground-based photometry (such 
as varying airmass and differing weather conditions). In- 
deed, it is clear that removal of systematics is a vital part 
of any transit search. The algorithm works without any 
previous knowledge of the particular systematics that af- 
fect the observations, and removes common trends in the 
data. The difference of each data point from the mean of 
the lightcurve is found, and the best linear fit determined. 
The slope is then found and the effect subtracted out, as- 
suming the same effect is present in many lightcurves. 

The algorithm does not increase the photometric preci- 
sion (as seen in Fig. [1]) except for the bright stars, whose 
scatter is dominated by systematics. For these the pho- 
tometric precision was increased from 0.01 magnitudes to 
0.003. From experimentati on with the transit s earch algo- 
rithms, we found that the iTamuz et al.l (j2005[ ) algorithm 
reduced the number of detected 'false transits' in the data 
by a factor of five, while leaving transit recoverability of 
true artificially-induced transits unaffected. 



5. PLANETARY TRANSIT DETECTION ALGORITHMS 

The transit search was carried out using two sepa- 
rate transit detection algorithms in tande m, namely the 
Box fitting Least-Squares (B LS) method of Kovacs e t al.l 
(2002) and the WS method of lWeldrake k Sa ckctt {200^. 
The WS and BLS codes themselve s are described in 
detail by IWeldrake fc SackettI (|2005f) and [Kovacs et"al] 
(|2002f ). respectively, to which we direct the reader for 
more detailed information. In general, the matched fil- 
ter approach for transit detection was first suggested by 
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iJenkins. Dovle fc CullersI (|1996( ) and has been described 
as the method of choice for t ransit searches in the gen- 
eral hterature (jTinglevI l2003al lbl) . Several recent experi- 
ments have adopted t he method w i th so me success (ie, 
iGilliland et all (|2000h : iBruntt et all (l2003h ). 

The WS version of the matched filter method assumes 
a simple square-well approximation for the transit shape 
(a justified assumption when searching for signals near the 
noise level) and involves comparison of the data to a large 
database of suitable transit models (~10^) with varying 
period, depth, duration and transit time. Each model 
is compared to the data via a stat istical test (a cross- 
correlation-function in the case of IWeldrake fc SackettI 
l2005f ) and a high significance output determines if a sig- 
nal has been detected and at which period. By intro- 
ducing detection criteria the significanc e of the detec- 
tion 'S jPniod : Tshift ) ' 1 as described in Wel drake fc Sacketj 
(|2005[ ). and the number of output data points above a 
pre-determined threshold (Np), the number of false detec- 
tions can be minimized while keeping the real recoverabil- 
ity high. 

The BLS method has also been utilized by many cur- 
rent transit surveys. The algorithm searches the data for 
the specific shape of a transit via the least-squares fit- 
ting of step-functions to the phase-wrapped data. This 
is repeated for a range of trial periods, fractional transit 
lengths, transit depths and transit epochs. The deviation 
of the fit to the data is calculated and the best-fitting 
combination of these parameters is flagged with a detec- 
tion significance. This significance (SDE) is somewhat 
dependent on dataset properties but has been generally 
set at in previou s apphcations (jKovacs et al.l 12002 ; 
iMocheiska et al.ll2006[ ). 

Both algorithms have similar transit recoverability lev- 
els as determined via separate Monte Carlo testing on our 
data. The WS code was implemented with strict transit 
detection criteria, S'(Pmod, Tshift)>ll and Np>10, dataset 
depen dent values as described in IWeldrake fc Sacketl 
(|2005[ ) that incorporate information about the real noise 
in the data and the window function. The BLS code, 
however, produces more accurate periodicity information 
for candidates, leading to easier identification when phase- 
wrapping and searching the candidates by eye. 

The degree to which the WS and BLS codes discriminate 
between transit-bearing artificial lightcurves and the total 
real dataset can be seen in Fig. |H The top histogram 
shows the BLS output SDE (cr) significance that results 
from running the algorithm on all dataset stars (open his- 
togram) and on only those stars for which an artificial 
transit has been added via our Monte Carlo simulations 
(shaded histogram). The difference in the mean signifi- 
cance between these two populations is only marginal with 
a large degree of overlap between the two distributions. 
Hence if BLS alone was used on the dataset, most, if not 
all of the stars, would need to be searched with another 
method to identify real transit features. The apparent gap 
in both distributions at 5a is statistical, being a product 
of our dataset window function and the chosen bin-size. 

For comparison, the bottom panel of Fig. [1] shows the 
WS output significance (S'(Pmod, Tshift)), both for the to- 
tal dataset (open histogram) and those with artificially- 
induced transits (shaded histogram). The difference is 



much more pronounced, with far less overlap, the transit- 
bearing stars having significantly higher output signifi- 
cance. By fiagging candidates above a certain detection 
significance (marked as a solid line) far fewer stars must 
be further analysed to determine real transit candidates. 
This threshold was set at 5'(Pmod, 7'shift)>ll for this data, 
as determined from tests to minimise false detection rates 
while keeping the recovery of transit features high. Both 
algorithms were therefore used in tandem, with WS used 
as a first-pass filter to produce a short list of candidates 
that were then further investigated with the BLS. 

All stars (both cluster and field) with suitable photomet- 
ric precision for the transit search were analysed with WS, 
producing an output candidate list. These candidates were 
then subjected to analysis with BLS to produce more accu- 
rate period information, hence making the transits easier 
to see when examined by eye. No further cuts were made 
on these candidates. They were then visually examined 
for transit features by eye, both at the peak periodicity 
and at integer aliases thereof. 

6. MONTE CARLO SIMULATIONS 

Extensive simulations with modeled transits of appro- 
priate depth and duration were carried out with the WS 
and BLS algorithms to determine the transit detection ef- 
ficiencies and hence expected transit recoverability in the 
dataset. From these recoverability results, the expected 
number of detected planets can be determined and in- 
formation can be gained on false-detection probabilities, 
which should be minimised as much as possible. 

A sample of 15,307 fake transit lightcurves were pro- 
duced using_j£tuai dataset lightcurves (after application 
of the iTamuz et al.l (|2005[ ) systematics removal package) 
for each of three photometric uncertainty bins: stars 
with precision <0.01 mag rms, 0.01<mag<0.02 rms and 
0.02<mag<0.04 rms. Over all three bins we have a total 
of 183,684 fake transits to test the detection algorithms. 

A box-shaped transit is superimposed on the lightcurve 
with a depth and duration consistent with an orbiting 
planet of various radii, assuming the star to be a clus- 
ter member with a previously determined radius (seen 
in Table. 1). Transit models incorporating stellar limb- 
darkening were not produced as this effect would not affect 
transit visibility and/or recoverability in our data. 

For the first bin (foreground stars) the stellar radius was 
assumed to I.OR0. For the two cluster bins, the stellar 
radius was determined as the weighted mean stellar ra- 
dius found from the cluster theoretical isochrones. These 
radii are 1.23Ro and O.OOR© respectively. For each depth, 
duration and orbital period bin, each model assumes a 
different epoch for the transit, spread randomly over the 
length of the observing run. By superimposing transits on 
randomly chosen dataset lightcurves, the models have the 
temporal resolution and photometric accuracy inherent in 
the real data. For the first bin (foreground stars), planets 
with radius l.ORjup and l.SRjup were simulated, whereas 
for the two cluster bins transits were produced that mimic 
the signal of planets with radius of only l.SRjup. We are 
insensitive to smaller cluster planets. 

The highest precision bin is appropriate only for fore- 
ground disk stars (as the other stars of this brightness 
would be cluster giants), whereas the other two bins sam- 



Planet Frequencies in w Centauri 



7 



pie the main sequence of the cluster. The artificial tran- 
sits were generated over 131 orbital periods incremented 
in 0.05 day steps in the range 0.52^7.10d and 117 tran- 
sit phases incremented by 0.2 days spanning the full MJD 
range of the dataset. The possibility of subday periods 
and a cla ss of 'Ultra Short Period' planets has been sug- 
gested by IS^ueTnjdOOi. An upper orbital period limit 
of seven days was used for the simulations, as this is the 
maximum period for which three transits may be visible 
over the observing run. Analysis of these artificial transits 
with both WS and BLS allows the sensitivity to planets 
of various radii to be determined, including the effects of 
the observational window function. 

All artificial transits were first analysed with WS, which 
produced a list of high significance (^(Pmod, 'n5hift)>ll) 
preliminary candidates. These candidates were then anal- 
ysed with BLS for more accurate period determination. 
Those recovered transits all had BLS-determined periods 
within ibO.ld of the real injected period. 

Examples of artificial transits that were detected at the 
low end of the detection criteria (SDE ~6.0(T significance 
with BLS, 5'(Pmod, Tshift)>ll with WS) can be seen in 
Figini indicating the visibility of transits in the dataset 
when phase-wrapped to the peak detected period. The 
panels are arranged in order of increasing photometric un- 
certainty, with typical orbital periods of VHJ (1.5d, left 
panels) and HJ (3.3d, right panels). 

The top two panels of Fig. [S] show a transit (at $ =0 
and 1.0) for a l.ORjup planet transiting a star of 1% pho- 
tometry, typical for the transit visibility for the brighter 
foreground stars. Moving downward, the two pairs of pan- 
els show a transit of a 1.5Rjup planet and a star of 2% pho- 
tometry, and a transit of l.SRjup against a 4% star, which 
defines the bottom of our search regime. The bottom two 
sets of panels are appropriate for stars on the cluster main 
sequence. 

Fig. [6] shows the resulting transit recoverability (as 
a fraction of the total number of modeled transits per 
period bin) for those transits for stars with <0.01 mag 
rms (foreground stars, top left) and 0.01<mag<0.02 rms 
(top right). These transits are those that passed both 
the S'(Pmod, 'n5hift)>ll and Np >10 detection criteria of 
the WS, and had periods within O.ld of the true injected 
transit period (although this in itself is not a criteria that 
a candidate must pass). These criteria were identical to 
those applied to the real data. For the brightest stars it 
can be seen that the recoverability is good to Porb=2d for 
those planets with radii >1.0Rjup. A large drop in recov- 
erability can be seen for periods approaching one day for 
all planetary radii, caused by terrestrial effects. 

As the photometric precision of the data becomes worse, 
only transits corresponding to larger planets can be de- 
tected. For stars with 0.01^0.04 mag rms, the recover- 
ability of planets with R < 1.5Rjup becomes increasingly 
truncated. By considering this, as well as the expected 
transit depths as a function of V magnitude in the cluster 
(see Fig. [Ij , the search can only provide meaningful statis- 
tics on planets with a radius of 1.5Rjup- The chosen upper 
limit is arbitrary, but is consistent with the upper limit of 
currently known transiting planets, par ticularly the dis- 
coveries of HAT-Pl- b and WASP-lb (Ba kos et al.ll2007at 
ICameron et al.|[2007( ) that are approaching this radius. 



The bottom left panel of Fig. [5] displays the l.SRjup 
transit recoverability for stars with photometric precision 
between 2— >4% (0.02— >0. 04 mag). These stars correspond 
to those with V=18.5^19.5 on the cluster main sequence. 
This defines the lower magnitude limit of the search, as 
the transit depth becomes smaller than the photometric 
scatter below this point. The recoverability for these stars 
is good for Porb <1.5d. The bottom right panel of Fig. 
[B] displays the weighted mean transit recoverability for all 
simulations, taking into account the relative numbers of 
stars that lie in each photometric bin. These numbers 
were used as part of the determination of the expected 
number of planets detectable in the cluster. 

Fig. [7]shows the transit recoverability for 1.5Rjup plan- 
ets as a function of V magnitude on the cluster main se- 
quence, assuming that they have short periods. The left 
panel depicts the recoverability of VHJ planets (with an 
assumed orbital period of 1.5d); the right panel shows the 
same for HJ planets (assumed period of 3.3d). The rapidly 
increasing stellar radii past the MSTO defines the bright 
limit to the recoverability, whereas the increasing photo- 
metric uncertainty (Fig. [1]) drives the faint limit. It can 
be seen that for VHJ, the recoverability is good for stars 
18.0<V<19.0. For HJ planets with R < 1.5Rjup, the re- 
coverability is low, as this is the lower limit to the plane- 
tary radius for which we can search at such faint magni- 
tudes (Fig. [1]). 

7. SENSITIVITY TO SHORT-PERIOD PLANETS 

With the results of the Monte Carlo simulations, we 
now determine our sensitivity to planets in uj Cen consid- 
ering a number of factors, including the total number of 
stars monitored with the appropriate photometric preci- 
sion, the transit probability and transit recoverability (as 
determined via the Monte Carlo simulations) as a func- 
tion of planetary period, and the assumed distribution of 
planets across orbital period. This sensitivity allows us 
to assess the significance of the results from our planet 
search, which we conduct for short-period planets with or- 
bital periods less than seven days, and thus encompassing 
both Very Hot Jupiter (VHJ) and Hot Jupiter (HJ) exo- 
planets. The essential question to ask is, "For how many 
stars would we have detected a transit signal if they were 
hosts to short-period planets?" In this section we answer 
that question, and tabulate the result under different as- 
sumptions in Table 2. 

Since both the probability that an exoplanet will tran- 
sit from the viewer's perspective and the probability that 
it will be detected with our transit algorithm are de- 
pendent on the planetary orbital period, an assumption 
must be made about how VHJ and HJ planets are dis- 
tributed across period. We consider four different sce- 
narios, VHJ planets distributed equally over the range 
1.0<P<3.0d (VHJ); HJ planets distributed equally over 
either 3.0<P<5.0d (HJ5) or over 3.0<P<7.0d (HJ7); and 
short-period planets distributed equally over 1.0<P<5.0d 
periods (VHJ+HJ). For each hypothesis, the 31,000 stars 
in to Cen for which we have suitable precision to detect 
transiting planets are uniformly apportioned to the rel- 
evant period bins, and given in column 5 of Table 2 as 
N 

Not all of these planets will transit their host stars, so 
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that the transit probabihty must be factored in. This 
probabihty, Ptrans,i7 depends on the planetary orbital pe- 
riod and radius of the host star, via: 



) 



(3) 



where i?* is the stellar radius and a is the planetary semi- 
major axis. 

The radius of the host star (i?*) varies with V magni- 
tude in the cluster, hence a weighted mean radius for the 
transit search stars was used in the calculation. By con- 
sidering the stellar radius from 17.5<V<19.5 (Table 1) in 
bins of half a magnitude, as well as the relative fraction of 
the total sample present in each bin, the weighted mean 
stellar radius was calculated to be O.96R0. The transit 
probability was then found for each planetary orbital pe- 
riod Pi, and is displayed in column 3 of Table 2. 

The final factor to consider is the probability that our 
algorithm will detect or recover the transit. Determined 
from the Monte Carlo simulations, this recoverability, 
R.trans,i, incorporates the window function and the true 
noise characteristics of the data. It is calculated for each 
planetary radius and the middle of each period bin from 
the weighted mean transit recoverability (bottom right 
panel of Fig. [S]). These values are given in the fourth 
column of Table 2 assuming a planetary radius of 1.5Rj. 

The total number of equivalent stars, N^^pj-obc, probed 
at full sensitivity for planets distributed according to one 
of our four hypothesis is then given by the simple sum 
of the numbers recorded in column 6 of Table 2, ie, 

N*, probe =I]jN*,inon,iXPtrans,iXRtrans,i- Fig. [5] showS this 

proxy for sensitivity, Nprobc,!, as a function of orbital pe- 
riod for each of the four hypotheses and assuming a plan- 
etary radius of l.SRjup in each case. Smaller planets will 
escape detection in our uj Cen experiment. 

8. SEARCH RESULTS 

Using the same search methods as in the Monte Carlo 
simulations, all 45,406 (31,000 cluster) stars with suitable 
photometric accuracy (<0.04 mag) in our lo Cen dataset 
were first analysed with the WS algorithm, and a candi- 
date list produced from those lightcurves that displayed a 
<S'(^'mod, Tshift)>ll output significance and a Np>10 crite- 
ria from the WS algorithm (see Fig. 2]). 

After filtering for detections with common inherent fea- 
tures (ie, those detections that occur at the same times 
and with the same periods, within O.OOld), 138 candi- 
dates remained. These candidates were analysed with BLS 
and their peak periodicities determined. They were then 
phase-wrapped to this period (and integer aliases) and vi- 
sually examined for transit features. 

The nature of the periodicity was immediately seen from 
the phase wraps. The vast majority of the candidates 
had transit-like features that could be attributed to bad 
columns in the CCDs, to non-perfect PSF fitting to the 
star due to blended companions, or association with stars 
close to the dataset saturation limit. Some candidates 
were composed of single outlying data points and others 
of random groupings of data points close to integer days 
which had no apparent pattern. Candidates with such fea- 
tures were classified as false positives. 



After this culling, none of the candidates (either in the 
foreground or the cluster) displayed transits of suitable 
depth and duration that could be attributed to an orbit- 
ing giant planet in our radius regime. The list did, how- 
ever, include 8 eclipsing binaries and 2 5 Scuti stars that 
were pr eviously identified and published in I Weldrake et all 
(|2007bD . 

The lack of planetary detections cannot be attributed to 
our algorithms, as evidenced by the Monte Carlo simula- 
tions and a control field in the Lupus Galactic Plane that 
we observed for 53 nights during 2005 and 2006 with the 
same instrument and detector. The Lupus data were an- 
alyzed with the same methods as used for lo Cen. Several 
transit candidates were found in Lupus, including one can- 
didate whic h has excellent prospe cts for being a transiting 
HJ planet (j Weldrake et al.ll2007f ). No similar candidates 
were seen in lo Cen. 

The net result of our search for transiting planets in 
u Cen is thus: no planetary candidates, 10 variable stars 
(previously detected with our variable star search algo- 
rithm), and 128 false positives. In the next subsection, 
we examine the implications of our null detection of large, 
short-period planets in the cluster. 

8.1. Upper Limits to Planetary Frequency in lo Cen 

An upper limit to the occurrence frequency F of large- 
radius VHJ and HJ planets in a; Cen can be determined 
using the results of our survey. Using Poisson statis- 
tics to analyze the significance of our null result under 
each of our four distribution hypotheses, we place an up- 
per limit on the fraction of stars in w Cen that are or- 
bited by short-period planets. For example, VHJ planets 
with 1.0 < P < 3.0d must have a frequency of occurence 
^VHJ < 1/1040 in order for none to be detected 95% of 
the time in our sample of N^.piobo =3100 "equivalent stars" 
probed at full sensitivity. 

Similarly, at 95% confidence, the frequency of HJ plan- 
ets in w Cen is i^HJs < 1/150 or Fhj7 < 1/93, depending on 
whether they have periods P that are equally distributed 
over 3 — 5 days or 3 — 7 days, respectively. If there is 
no strong difference in occurrence frequency between VHJ 
and HJ, so that we may assume that the combined pop- 
ulation is equally distributed over Id < P <5d, the up- 
per limit on the frequency of such short-period planets is 
^VHJ+HJ < 1/600 (95%CL). These upper limits are tabu- 
lated in Table 2. 

8.2. Metallicity considerations 

In the Solar Neighborhood, evidence suggests that plan- 
etary frequ encies are influenced by the metallicity of the 
host star (iGonzalez 1997;_ Laughlinl I2OOOI : ISantos et al.l 
[200lt [Fischer fc Valentill2005fl . as also appears t o be the 
case in the globular cluster 47 Tucanae (Weldr ake et al.l 
l2005f ). Based on statistics gathered in radial velocity 
searches, iMarcv et al.l (|200l) conclude that the probabil- 
ity, Ppianet, that a given star will host a planet is related 
to metalhcity via Ppia„ct =0.03x 102 °P<=/h]^ where [Fe/H] 
is the stellar metallicity and the constant of proportion- 
ality is assumed to be the same for planets of all orbital 
periods. At very low metallicities, the relation is poorly 
constrained observationally, and thus highly uncertain. 



Planet Frequencies in w Centauri 



9 



Of the 31,000 monitored cluster stars in lu Cen, 80% 
are assumed to be members of the main metal-poor pop- 
ulation, w ith co rresponding [Fc/H]=— 1.7, taken from the 
results of iNorris (2004). Similarly, 15% and 5% are as- 
sumed to have [Fe/H]=— 1.2 and —0.6 respectively. The 
weighted mean metallicity for the sampled stars is hence 
[Fe/H]= -1.57. 

If the above relation between Ppianct and stellar metal- 
licity also holds true for globular clusters, then the proba- 
bility that planets will be found around stars in lu Cen is 
reduced, by an amount that depends on the extrapolation 
of the relation into uncertain territory. A conservative es- 
timate, that we will assume in the next section, is that 
w Cen planetary occurrence frequency is sup pressed by at 
least t he ^60% reduction in F estimated bv iMarcv et al.l 
()2005[ ) for exoplanets around stars of [Fe/H]=— 0.5 com- 
pared to those with nearly solar metallicities. In the next 
section, we compare our results with rates and frequency 
limits of short-period planets in non-cluster Galactic envi- 
ronments. 

9. DISCUSSION AND IMPLICATIONS 

We have placed upper limits on the frequency with 
which Lo Cen stars host short-period (VHJ and HJ) plan- 
ets. How do these limits compare to the frequency of VHJ 
a nd HJ planet s detec ted in other Galactic fields? 

iGould et af] (|2006t) analyzed the OGLE HI transit sur- 
veys in Galactic disk and bulge fields to derive occur- 
rence rates for VHJ and HJ planets with radii satisfying 
1.0 < Rjup < 1.25. These rates, according to the au- 
thors, are statistically indistinguishable from those derived 
from radial velocity results on bright nearby stars. No 
larger, short-period planets were detected by t he OGLE HI 
survey s despite increased sensitivity, leading [Gould et all 
to compute an upper limit F for these larger worlds. 

For Very Hot Jupiter (VHJ) planets unif o rmly dis- 
tributed over Id < P < 3d, the iGould et all (|2006f ) re- 
sult is F =l/llO{ltl\l) at a 90% confidence level. Their 
95% confidence rate is thus somewhat larger than our up- 
per limit of FvHj < 1/1040 for somewhat larger planets 
in w Cen, possibly indicating suppression in the cluster. 
However, within the 90% confidence envelope, the two are 
consistent regardless of whether a correction for metallic- 
ity is applied or not. 

For Hot Jupiter (HJ) planets uniformly distributed over 
3d < P < 5d. IGould et"an (|2006D find F = 1/320(1^^^^ 
at a 90% confidence for smaller transiting planets in 
OGLE III fields. For larger planets in cu Cen, our upper 
limit is i%j5 < 1/150 using the same assumption about 
period distribution. Since the upper limit on HJ planets 
in this cluster environment is larger than the occurrence 
rate inferred for OGLE fields, even without a metallicity 
correction, no conclusion can be drawn. 

Perhaps the most direct comparison that can be made 
of the two studies is the upper limits on the frequency of 
i? = 1.5 Rj planets assuming uniform distribution over the 
orbital period range Id < P <5d. Our result for uj Cen is 
-PVhj+hj < 1/600. The corresponding value for OGLE HI 
fields can be inferred from Table 6 of IGould et all (|2006[ ) 
to be a very comparab le upper limi t of F < 1/640. 

iFregeau et al.l ()2006f ) (and earlier. ISigurdssonlll992 D car- 
ried out a study of dynamical interactions of planet- 



bearing stars with other stars in dense stellar environments 
via N-body simulations. They conclude that planet sur- 
vivability is greater in dense systems than had been pre- 
viously predicted, and presented a characteristic timescale 
for planet survivability for various star-planet-star mass 
ratios, systemic stellar densities and v elocity dispersions 
(see Fig. 6 and Eq 22 in lFregeau et al.| [2006). 

Lo Cen has a measured core stellar density of only 
1/10**^ the core density of 47 T uc, and a measured core 
velocity dispersion of 20km/s (|Mevlan fc Mavoij 119861 : 
Ivan de Ven et al1l2006f ). which falls to 14km/s at the clus- 
ter half-mass radius. Using these parameters, the expected 
lifetime for a planetary system to remain bound in the 
cl uster is >lx 10^° y ears. 

ISigurdss"onl (119921 ) has shown that the survivability is 
long (~10^yrs) for a planet of ~1AU semi-major axis in 
an environment typical of the core of 47 Tuc, and even 
longer for the lower density environment of u Cen. For 
the short period HJ planets, the probability of disruption 
due to stellar interaction is very low. The transiting planet 
search presented here is most sensitive to the outer halo 
of the cluster, where stellar densities are far lower than 
in the core. Hence 1x10^" years can be taken as a firm 
lower limit to the timescales of planet survivability in the 
cluster. 

The conclusion is that if the planets formed in the first 
place, then they are expected to survive in lu Cen through 
its current and remain detectable. The low upper limits 
on the frequency of short-period planets in the cluster are 
consistent with its low metallicity inhibiting planet forma- 
tion from the outset. 

Our null result does not rule out the existence of large- 
radius planets in longer period orbits in the cluster. The 
confirmatio n of a p lanetary mass object in the globular 
cluster M4 (jSigurd sson et al. 2003) provides the first di- 
rect evidence of planetary formation in a very metal poor 
environment. The paucity of transiting short-period large 
planets in globular clusters is due to a process other than 
stellar dynamics, perhaps a dependency of planetary mi- 
gration on stellar metallicity. 

10. SUMMARY AND CONCLUSIONS 

We have presented the results of a wide-field, deep 
photometric search for transiting short-perod planets in 
the globular cluster uj Centauri, a region previously un- 
sampled for planetary transits. The cluster was observed 
with a 52'x52' field of view for 25 contiguous nights with 
the ANU 40-inch telescope at Siding Spring Observatory. 
From application of difference imaging analysis, a total of 
109,726 time series were produced across the field, each 
being composed of 787 independent data points. 

A total database of 45,406 stars have photometric accu- 
racy suitable for the search (<0.045 mag scatter down to 
V=19.5), including 31,000 cluster stars extending 2.5 mag- 
nitudes down the main sequence. All of these were sub- 
jected to a rigorous (and vigorous) search for transit-like 
events; none were detected after variable stars and clear 
false-positives were removed. Simulations have shown that 
if large Hot Jupiters (HJs) formed in the cluster then dy- 
namically speaking they would survive to be detectable in 
our data. 

Extensive Monte Carlo simulations via injection of tran- 
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sit signals into actual light curves were used to determine 
the sensitivity of the survey to R<1.5Rjup planets over a 
range of orbital periods. Coupled with our null result, wc 
are thus able to place strict, statistically significant upper 
limits on the occurrence frequency F of large (R — 1.5R), 
short-period planets in u; Centauri. 

We determine a limit of Fynj < 1/1040 at 95% con- 
fidence for Very Hot Jupiter (VHJ) planets with periods 
distributed uniformly over Id < P < 3d. Thi s upper limit 
for the cluster is less than that determined bv lGould et al.l 
(|2006l ) for smaller (l.SRjup < R < 1.5Rjup) planets with 
the same period distribution in the Galactic fields sur- 
veyed by OGLE III. The two results are consistent at the 
90% confidence level, and more understandable if the low 
metallicity of lo Gen suppresses planet formation or plan- 
etary migration. 

Under the assumption that there is no difference in oc- 
currence frequency for VHJ and HJ across the orbital 
period range Id < P <5d, we derive an upper limit of 
^VHJ+HJ < 1/600 in uj Gen. The corresponding result in 
the Galactic OGLE III fields for comparably-sized planets 
is an upper limit of F < 1/640. Both results are quoted at 
95% confidence. Our results are less dependent on model 
assumptions about the distance to the target population 
since the vast majority of stars in our fields are members 
of, and thus at the distance of, ui Gen. 

It is noteworthy that despite the fact that the OGLE HI 
campaigns monitored considerably more stars in better 



median seeing with more frames per field, our to Gen study 
produces a comparable upper limit on the frequency of 
large, short-period planets. While part of the reason may 
lie in the longer exposure times, denser sampling, and 
the use of different cleaning and detection algorithms in 
our survey, a large part of the difference is due to the 
small fraction of the more distant and obscured Galactic 
bulge stars (which constituted about a third of the total 
OGLE III sample) that can be meaningfully probed for 
transiting planets. 

This null result for VHJ and HJ planet s in ui Gen, cou- 
pled with the null result of 47 Tucanae (jWeldrake et al.l 
l2005f ) strengthens the evidence for the dominance of sys- 
tem metallicity over stellar interactions in determining 
short period planetary frequencies in globular clusters. 
At longer orbital periods stellar encounters may play a 
role in determining planetary frequencies. This is a re- 
sult aligned with current work on the metallicity trend 
of planet-bearing host stars in the Solar Neighborhood 
and N-body simulations of planets in dense environments. 
Such a metallicity dependence is one of the main predic- 
tions of the core accretion model of planet formation. 

The authors would like to thank Omer Tamuz for discus- 
sions on the removal of dataset systematics, Laura Stan- 
ford for help with the production of the theoretical cluster 
isochrones, and the referee for a very thorough and helpful 
review process. 
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Table 1 



V 


Mass(MQ) 


T(K) 


L(Lq) 


log{g) 


Mv 


R(Rq) 


N, 


17.0 


0.7636 


5500 


4.08 


3.6212 


+3.4 


2.24 


1300 


17.25 


0.7613 


5700 


3.21 


3.7848 


+3.65 


1.85 


1700 


17.5 


0.7579 


6000 


2.64 


3.9375 


+3.9 


1.55 


2200 


17.75 


0.7487 


6100 


1.89 


4.1312 


+4.15 


1.23 


2800 


18.0 


0.7426 


6100 


1.62 


4.2003 


+4.4 


1.13 


3200 


18.25 


0.7272 


6100 


1.18 


4.3167 


+4.65 


0.98 


3060 


18.5 


0.7186 


6000 


1.01 


4.3653 


+4.9 


0.92 


2900 


18.75 


0.6993 


5900 


0.76 


4.4516 


+5.15 


0.82 


2700 


19.0 


0.6897 


5900 


0.65 


4.4917 


+5.4 


0.78 


2500 


19.25 


0.6687 


5800 


0.49 


4.5629 


+5.65 


0.71 


2200 


19.5 


0.6575 


5700 


0.43 


4.5945 


+5.9 


0.68 


2000 


20.0 


0.6092 


5400 


0.26 


4.6978 


+6.4 


0.58 


1800 



Not e. — The determined stellar parameters for uj Cen turnofT and main sequence stars from V=17.0 to V=20.0 as produced from lYl et al.l 
l|2003l ) theoretical isochrones. The tabulated values are the apparent V magnitude (using the best-fit distance modulus of 13.6 and a metallicity 
of — 1.7dex), the stellar mass in solar units, the stellar temperature and luminosity in solar units, the logarithm of the surface gravity, the 
corresponding absolute magnitude and the determined stellar radius along with the estimated number of cluster stars present in each magnitude 
bin. These radius considerations were used in determining transit recoverability and visibility. 



Table 2 



Pi fdavs) 




p 

trans, 1 


rj 

tr&ns ; 1 


N 

* mon ji 


*proDe,i 
















0.5 





~ 








l.O 


0.25 


0.20 


0.84 


7750 


1302 


1.5 


0.25 


0.18 


0.70 


7750 


977 


2.0 


0.25 


0.15 


0.43 


7750 


500 


2.5 


0.25 


0.14 


0.31 


7750 


336 












N. nrohp =3100 

^Vhj <1/1040 


3.0 


0.25 


0.12 


0.19 


7750 


177 


3.5 


0.25 


0.11 


0.16 


7750 


136 


4.0 


0.25 


0.10 


0.13 


7750 


101 


4.5 


0.25 


0.09 


0.07 


7750 


49 












probe — 


3.0 


0.125 


0.12 


0.19 


3875 


88.4 


3.5 


0.125 


0.11 


0.16 


3875 


68.2 


4.0 


0.125 


0.10 


0.13 


3875 


50.4 


4.5 


0.125 


0.09 


0.07 


3875 


24.4 


5.0 


0.125 


0.09 


0.05 


3875 


17.4 


5.5 


0.125 


0.09 


0.05 


3875 


17.4 


6.0 


0.125 


0.08 


0.04 


3875 


12.4 


6.5 


0.125 


0.07 


0.02 


3875 


5.4 












N*, probe =280 

Fhj7 <1/93 
















0.5 













1.0 


0.125 


0.20 


0.84 


3875 


651 


1.5 


0.125 


0.18 


0.70 


3875 


489 


2.0 


0.125 


0.15 


0.43 


3875 


250 


2.5 


0.125 


0.14 


0.31 


3875 


168 


3.0 


0.125 


0.12 


0.19 


3875 


88.4 


3.5 


0.125 


0.11 


0.16 


3875 


68.2 


4.0 


0.125 


0.10 


0.13 


3875 


50.4 


4.5 


0.125 


0.09 


0.07 


3875 


24.4 



N,,probe =1800 

^Vhj+hj < 1/600 



Note. — Numerical parameters used lo calculate the total number of stars probed for the four assumptions described in the text, for the 
distribution of VHJ planets as a function of period. The tabulated parameters as described in section 7 arc the orbital period for bin i, Pi, 
the transit probability, Ptrans.i the detection algorithm transit rccovcrability, Rtrans,ii for that bin, the number of stars monitored in the bin, 
N*mon,ii and the resulting number of stars probed N, probe, i in Cen. Since the total number (31,000) of ui Cen stars monitored with sufficient 
photometric precision is estimated to two significant figures, N„probe,i similarly estimated. For each hypothesis, an upper limit (95% CL) F 
to the occurrence frequency of short-period planets is also given. 
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Fig. 1. — The measured photometric rms uncertainty after apphcation of systematics removal for 97,935 stars in the ui Cen field that were 
cross-identified with the CMD dataset, measured as the standard deviation as a function of V magnitude. Also plotted as the red points is 
the total noise contribution (star+sky+residual noise) calculated for each star independently from their respective star and sky fluxes. Also 
plotted is the location of the cluster MSTO (at V=17.5). Various percentage uncertainties are marked for clarity. It can be seen that the 
photometric uncertainty is 0.01 mag at the cluster turnoff, degrading to 0.04 mag at V=19.0. Also overplotted are the expected transit depths 
of a planet with R=1.5Rjup for comparison with the photometry. Stars as faint as V=19.5 can be searched for R=1.5Rjup planets with a 
transit depth equal to the photometric uncertainty. We are insensitive to planets with smaller radii. 

Fig. 2. — The color magnitude diagram for the upper main sequence of lo Cen. The boxes con tain the stars s earched for transits based 
on the depth of the expected signal compared to the photometric precision. Overplotted are three lYi et al.l II2003I) theoretical isochrones to 
describe the cluster population with [Fc/H]=— 1.7 (solid line), — 1.2 (dotted line) and —0.8 (dashed line), values as taken from[Norris (2003). 
Assuming the best-fitting distance modulus of 13.6, the radius and mass of the main sequence stars were found as a function of V, and the 
corresponding depths (in magnitude units) and durations (in hours) of an orbiting 1.5 Rjup planet were calculated and overplotted. The 
Monte Carlo simulations were carried out with these parameters. The duration was calculated assuming an orbital period of 3.3d, typical for 
HJs in the solar neighborhood. 

Fig. 3. — The extent of our WFI field of view, shown as ARA and ADec (in degrees) from the core of the cluster. The locations of the 
eight WFI CCDs are overplotted with the locations of the cluster core radius (inner ellipse), the cluster half-m ass radius (middle ellipse) and 
the position of 50% of the cluster tidal radius (outer ellipse). Cluster parameters are taken from IHarrisI II1996I) . 

Fig. 4. — A comparison of the degree to which the BLS (top) and WS (bottom) detection codes discriminate between lightcurves with 
artificially induced transits (dark shading) and no added transits (open histogram) superimposed on our actual observational light curves. 
The histograms of significances for the BLS code are similar whether transits are present or not, hence most, if not all, of the stars must 
be examined by eye at the peak detected period to identify the transits. The WS code, however, has more disparate histograms for the two 
populations, leading to far fewer false detections in the output candidate lists. The BLS code produces more accurate period determinations. 
Hence it was decided to use both algorithms in tandem: the WS to first filter the lightcurves for high significance candidates (with a threshold 
marked by the solid line and an Np criterion of >10 points) and the BLS peak detection for accurate period determination. 

Fig. 5. — Examples of the visibility of artificial transits producing output significances at the lowest end of detection significance (SDE ~6 cr 
with BLS, 11 ^(Pniodi ■^shift) with WS). Left panels show a transit for a planet with an orbital period (Porb) of 1.5d (typical of a VHJ) with the 
right panels having Po^b =3. 3d (HJ). From top to bottom are transits for a l.ORjup planet superimposed on a 1% photometry (foreground) 
star, a 1.5Rjup planet with a 2% star and a 1.5Rjup with a 4% star. The transits have been produced with parameters appropriate for the 
radius of the host star (Table. 1). 

Fig. 6. — The statistical transit recoverability for a range of planetary radii for stars with photometric precision <0.01mag (top left) 
and 0.01— >0. 02 mag (top right), 0.02^0.04 mag (bottom left) and weighted mean (bottom right) plotted as the fraction of artificial transits 
recovered as a function of orbital period in days. The top left panel corresponds to the brightest stars, hence Galactic disk contamination. 
The top right panel corresponds to stars on the upper 0.8 magnitudes of cluster main sequence. The drop in recoverability near Id is due to 
diurnal effects. These window functions indicate the necessity of transit surveys to contain many data points over a long observational period 
in order to be statistically sensitive to longer period planets. The transit recoverability for stars towards the faintest end of our search range 
(bottom left), typical of stars from V~18.2— >19.0. The bottom right panel shows the weighted mean transit recoverability from all Monte 
Carlo simulations across all stars, calculated accounting for the relative numbers of cluster main sequence stars in our search range. This 
recoverability was used in calculating the expected number of planets that should be detectable in our dataset. 

Fig. 7. — The recoverability fraction of planetary transits in our dataset for R=1.5Rjup planets as a function of V magnitude along the 
cluster main sequence. The limits to the recoverability are produced by the larger sizes of cluster stars above the cluster main sequence 
(V<17.5) and the increasing photometric uncertainty (V>19.5). 

Fig. 8. — The total number of stars probed for full sensitivity to transiting large radius planets in the cluster, distributed in orbital period 
for each of our four period distribution hypotheses. The solid line denotes the number of stars for VHJ planets distributed equally over 
1.0< P <3.0d, the short-dashed line shows the star numbers for HJ planets of 3.0< P <5.0d, with the long-dashed line for HJ planets of 
3.0< P <7.0d and finally the dotted Une for VHJ and HJ distributed equally over 1.0< P <5.0d. 
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